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Macroscopic fluctuation theory and first-passage properties of snrface diffnsion 

Baruch Meerson^’El and Arkady Vilenkin^ ’0 
^ Racah Institute of Physics, Hebrew University of Jerusalem, Jerusalem 9190f, Israel 

We investigate non-equilibrium fluctuations of a solid surface governed by the stochastic Mullins- 
Herring equation with conserved noise. This equation describes surface diffusion of adatoms accom¬ 
panied by their exchange between the surface and the bulk of the solid, when desorption of adatoms 
is negligible. Previous works dealt with dynamic scaling behavior of the fluctuating interface. Here 
we determine the probability that the interface first reaches a large given height at a specified time. 

We also And the optimal time history of the interface conditional on this non-equilibrium fluctuation. 

We obtain these results by developing a macroscopic fluctuation theory of surface diffusion. 

PACS numbers: 05.40.-a, 05.70.Np, 68.35.Ct 


The stochastic Mullins-Herring equation with con¬ 
served noise can be written in a rescaled form as 


dth = -V'^h + V (1) 

where h(x,t) is the interface height, and £ is a Gaussian 
noise which is delta-correlated in x and t and has zero 
mean and unit variance. 

The noiseless version of Em o was suggested by 
Mullins almost 60 years ago [l[. It describes the cap¬ 
illary flattening of a perturbed solid surface to its equi¬ 
librium shape, where the surface diffusion of adatoms is 
accompanied by adatom exchange between the surface 
and the bulk of the solid, while the adatom desorption 
is negligible Equation © was extensively studied 

in the context of dynamic scaling behavior of fluctuating 
interfaces lM3l- A closely related discrete model is the 
restricted solid-on-solid model with conserved noise m- 
In this model a pair of nearest-neighbor sites {i,i + I) 
is randomly selected on a (d — I)-dimensional substrate. 
One particle is moved from the site i to the site i + 1, or 
vice versa, with probability 1/2. The move, however, is 
only allowed if the inequality \h{i -|- l,t) — h(i,t)\ < 2 is 
satisfied after the move. 

The previous works 0, H, [HI focused on the inter¬ 
face width w{L,t) which, for the conserved equation like 
Eq. (P), is defined by w‘^{L,t) = {h?{v,t)), where L is 
the linear size of the system. Assuming an initially fiat 
interface, w{L, t) exhibits scaling behavior: 


r , 0 < f 

I 


( 2 ) 

( 3 ) 


where /3 = (2—d)/8, a = (2—d)/2, and d is the dimension 
of space [7Hll|. The dynamic exponent z = 4 is related 
to a and /3 via z = a//?. The most interesting case here 
is d = 1, when both the growth exponent /3 = 1/8, and 
the roughness exponent a = 1/2 are positive. That is, in 
one dimension the small intrinsic noise makes the inter¬ 
face roughen without any external driving, in spite of the 
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smoothing effect of the surface diffusion. The roughening 
disappears if one replaces in Eq. o the surface diffusion 
term —V^h by the more customary term V^h [si-fTol. [l^ . 

The interface width w{L,t) is a useful integral mea¬ 
sure of the fluctuating interface. Here we will address 
Eq. (P from a different angle, by focusing on a first- 
passage property of the interface [Q. We will evaluate 
the probability that the interface first reaches a (large) 
given height at time T, and find the optimal history of 
the interface conditional on this non-equilibrium event. 

To achieve this goal we will develop a macroscopic fluc¬ 
tuation theory (MET) for Eq. p. This is a WKB-like 
theory in the spirit of the (weak-noise limit of) Martin- 
Siggia-Rose field-theoretical formalism H- More re¬ 
cently, variants of this theory were developed for diffusive 
lattice gases: for their non-equilibrium steady states, see 
fisl for a review, and for non-stationary settings [H- 
ll8|. Using the MET formulation, we will first evalu¬ 
ate the probability to observe, at t = T, a specified 
height profile hT{x) when starting from a fiat profile at 
t = 0. The solution also gives the optimal height profile 
history leading to this hrix). Then we will deal with 
the first-passage problem. Here there are two different 
regimes: the equilibrium regime at T ^ L^, and the non¬ 
equilibrium regime at T In the equilibrium regime 

the first-passage probability can be found from a mini¬ 
mization of the free energy of the system under proper 
constraints, whereas the optimal activation history co¬ 
incides with the time-reversed relaxation history, as to 
be expected from the Onsager-Machlup symmetry M- 
In the non-equilibrium regime a full time-dependent solu¬ 
tion of the MET equations is needed for the evaluation of 
the first-passage probability, whereas the activation and 
time-reversed relaxation histories are entirely different. 
In the non-equilibrium regime we uncover dynamic scal¬ 
ing behavior of the first-passage probability and of the 
height profile corresponding to reaching the maximum 
height at t = T. Finally, we will show how to extend 
some of our results to the Mullins-Herring equation with 
non-conserved noise. 

Let us measure the coordinate x along the substrate, 
|a;| < T/2, and assume periodic boundaries. The MET 
equations follow from a saddle-point evaluation of the 
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action integral, obtained in a standard way from Eq. (HD 
[^ . They can be written as two Hamiltonian partial 
differential equations: 

dth = 6H/Sp= -dxid^h + dxp), (4) 

dtp = -SH/Sh = d*p, (5) 

with the Hamiltonian 

7? = y dxH, 'H = dxP + ^9xP^ ■ (6) 

With proper initial and boundary conditions, that we 
will discuss shortly, Eqs. (HJ and ([5]) describe the optimal 
interface height history h{x, t). Once they are solved, we 
can evaluate the probability V of the large deviation we 
are interested in: — In P ~ S', where 

pT pL/2 

S= dt dxipdth — H) 

Jo J-L/2 

1 /•^/2 

= - dt dx{dxpf (7) 

^ Jo J-L/2 

is the action. Let us briefly discuss some general prop¬ 
erties of this Hamiltonian flow. As follows from Eq. m, 
there are two invariant zero-energy manifolds. The mani¬ 
fold dxP = 0 corresponds to the noiseless Mullins-Herring 
equation dth = —d'^h. The second invariant zero-energy 
manifold, described by the equation 

dxP = -2dlh, ( 8 ) 

corresponds to thermal equilibrium. Using Eqs. o and 
dH), we obtain 

dth = dxh. ( 9 ) 


For finite T the system is out of equilibrium, and we 
need to solve Eqs. (SD and m explicitly. We can write 

OO 

h{x,t) = ^ a„(t) cos(fc„a;)-I-6„(t) sin(fc„x), (11) 

n—1 

OO 

p{x,t) = '^an{t)cos{knx) + I3nit)sm{knx), (12) 

n—1 


where = 27rn/L. Solving the ensuing ordinary differ¬ 
ential equations for the coefficients, we obtain 


a„(t) 


'A 

bn{t) 


B 

CXn (t) 


'A 

Pn{t) 


B 


sinh(fc^t) 
sinh(fc)(T) ’ 

sinh(fc((T) ’ 


and 


(13) 

(14) 


where A„ and are the Fourier coefficients of hx^x). 
Now we evaluate the double integral in Eq. © and obtain 


S = LY,KniL,T){Al + Bl), (15) 

n—1 


where 

= 2(1-%tTy 

Equations m and (fTBl) yield the probability of observing 
a given height profile hx(x) at an arbitrary time T . This 
is a simple but fully non-equilibrium result. 

For T Kn{L,T) ~ k^/2, and Eq. dTSl) reduces to 
the free energy from Eq. dloD, whereas the activation tra¬ 
jectory coincides with the time-reversed relaxation tra¬ 
jectory. For T the system is far from equilibrium. 

As a simple example, suppose that hx^x) only includes 
long wavelengths, comparable to L. Then 


That is, the equilibrium dynamics is described by the 
time-reversed Mullins-Herring equation. As a result, an 
activation trajectory at thermal equilibrium coincides 
with the time-reversed relaxation trajectory [l9j |. 

Before dealing with the first-passage problem, we will 
solve an auxiliary problem by specifying a height profile 
hx{x) at time t = T. For simplicity, we assume a flat 
interface at t = 0 and set h{x, t = 0) = 0. In the limit of 
T /$> L'^, the system will explore equilibrium fluctuations 
in order to reach hx{x). In this limit there is no need 
to find the activation trajectory: It suffices to evaluate 
the difference between the free energies of the final and 
initial states. Indeed, let us evaluate the action © on 
the equilibrium manifold ([5|), using Eq. (0: 


S = 2 J dt J dx{dxhY = 2 J dt J dxd^hdxh 


= 2 dt dxd‘l^hdxh = 


rL/2 


l-L/2 


dx [dxhx(x)f-.{IG) 


q _^ 

~ AT ^ kl ’ 

which diverges at T —>■ 0 implying a very small prob¬ 
ability of observing a large given height at short time. 
In their turn, the Fourier coefficients an and bn from 
Eq. (fT51) can be approximated as a„(t) ~ (t/T)A„ and 
hn{t) — {t/T)An, so that the optimal interface height 
history becomes simply h{x,t) — {t/T)h{x,T). This is 
very different from the time-reversed relaxation history 
at equilibrium. 

With Eqs. (fTT]) - (IT6l) at hand, we can now evaluate the 
probability that the interface first reaches a (large) given 
height, say M, at time T. Without losing generality, 
we assume that this happens at a; = 0. We start from 
minimization of the action (1151) under constraint hx{x = 
0) = M, without demanding that this height was not 
reached at earlier times 0 < t < T. In view of Eq. (HD, 
the constraint hx{x = 0) = M becomes 

OO 

J2An=M. 

n—1 


The final expression is the free energy cost of the height 
profile hx{x), as to be expected. 


( 17 ) 
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Let us introduce constrained action, 


5a = ^ [LK4L,T){Al + Bl) - AA„] , 

n—1 

where A is a Lagrange multiplier. Varying S\ with re¬ 
spect to An and Bn, we obtain 

OO 

^ [2LKn{L,T){AndAn + BnSBn) - AM^] = 0, 

n—1 


which yields 


A 


n 


A 

2LKn{L,Ty 


Bn = Q, n = 1, 2,.... 


Plugging these An into Eq. (IT71) . we find 

2LM _ LM 
~ Y.n=iKn\L,T)~Q{i^y 


where 


Q{L,T) 


E 


1 — e 


-2ktT 


1.2 


(18) 


(19) 


Now we use Eq. (jl3l) and find the time-dependent Eourier 
coefficients a„ and bn'. 


an{t) 


M sinh(fc((t) 

2Kn{T)Q{L,T) sinh(fc4T)’ 


6„(t)=0. (20) 


Therefore, the optimal height profile history is 


h{x, t) 


M ^ 1 — e sinh(fc((t) 

kl sinh(fc4r) 


COs(knX). 


( 21 ) 

Importantly, the value h = Af is achieved at t = T for 
the first time. Therefore, Eq. m solves the first-passage 
problem that we are after. Figure [T] shows the optimal 
interface height histories h(x, t) for large (the top panel) 
and small (the bottom panel) values of T. As one can 
see, these histories are very different. This includes the 
height profiles at t = T, 


h{x, T) 


M ^ 1 - 
Q{L,T)^^ kl 


COs(knX). 


( 22 ) 


which, surprisingly, develop a corner singularity at x = 
0 where the value h = M is first reached. Note also 
that at large T (the top panel) most of the dynamics 
happens at times close to T: the system has enough time 
to thermalize before producing a large deviation. 

The first-passage action can be found from Eq. (USD: 


-laV~S = 


LM"^ 

2Q{wry 


(23) 




X 


FIG. 1; (Color online) The optimal interface height history 
conditional on reaching h — M at t — T. The initial state 
is h(x,t = 0) = 0. Shown is the rescaled height h(x,t)/M 
versus the rescaled coordinate 2'itxIL. Top panel; T = 5 and 
f = 3 (dotted line), 4.99 (dashed line), and f = 5 (solid line). 
Bottom panel: T = 0.01 and t=0.003 (dotted line), 0.009 
(dashed line), and t = 0.01 (solid line). Time is rescaled by 

(Li2^y. 


The resulting probability V{M) is Gaussian, with the 
standard deviation a(L,T) = Q{L, T)/L. The asymp¬ 
totics of the function Q(L,T) are the following: 


Q(L,T)^ 


L^/24, 

r(3/4)ri/4L 

2374^ 


T » L^, (24) 

T^L^, (25) 


where r(...) is the gamma-function. The T ^ 
asymptotic is obtained by neglecting compared to 

1; the remaining infinite sum is equal to L^/24. To obtain 
the T <C asymptotic, we replaced the infinite sum in 
Eq. dT51) by the integral which yields Eq. (1251) . The stan¬ 
dard deviation ol'P{M) has the following asymptotics: 


<7(L,T)~ 


2,1/2 

[r(3/4)]4/2T4/8 

23/87rl/2 


T » L4, (26) 

T <C (27) 


At long times a is independent of T, and its scaling with 
L is described by the roughness exponent a = 1/2. At 
short times a does not depend on L, and its scaling with 
T exhibits the growth exponent [3 = 1/8. 
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FIG. 2: (Color online.) The function ‘I’(^) = h{x,T)/M from 
Eq. dSO]). 


Now let us examine the optimal height profiles at t = T 
in the long- and short-time limits. In the long-time limit 


h{x, T 


oo 


6M ^ 

n—1 


COs{knX) 


= M 


ii) 


2 61x1 


-hi 


(28) 


This shape is very close to that observed at t = T on 
the top panel of Fig. [T] In the short-time limit we can 
replace the infinite sum in Eq. (l22l) by an integral. After 
a change of variable, 


(29) 

where ^ = xI{2TY'‘^. The profile h{x,T) is self-similar, 
as X and T enter only through The integral in Eq. (1^ 
can be expressed via the generalized hypergeometric 
function [2l|, and we obtain h{x,T <g; 1) = M$(^), 

V 4’ 4’ 2 ’ 4’ 256 J 


2r(!) 


.(30) 


The graph of 4>(^) is shown in Eig. [21 The “optimal” 
interface shape is oscillatory; it is very close to the one 
observed at t = T on the bottom panel of Eig. |TJ 

The long-time optimal height profile (1251) can also 
be obtained from a minimization of the equilibrium 
free energy m, constrained by the conservation law 

L /2 

f_j^/^hTix) dx = 0. Indeed, introducing a Lagrange 
multiplier A, we arrive at a simple minimization prob¬ 
lem for the Lagrangian 


L{hT,h'j,) = {h'j,f - Khr, (31) 


where the primes denotes the x-derivative. The Euler- 
Lagrange equation is 2h^ -I- A = 0. Its general solution 
is a parabola: 

^^(x) = Cl-I-C 2 X — (l/4)Ax^. (32) 


As hriO) = M, we obtain Ci = M. Importantly, the 
solution cannot be smooth, because one cannot satisfy 
three additional conditions - the conservation law and 
the periodicity conditions hT{—L/2) = hT[L/2) and 
h'rp{—L/2) = h'rp{L/2) - with only two arbitrary con¬ 
stants C 2 and A. The way out is to allow a discontinuity 
of dh/dx at X = 0, so that the maximum at x = 0 is 
non-analytic. In this case the coefficient C 2 does not 
vanish, and has the same magnitude but opposite signs 
at X < 0 and x > 0. This immediately leads to the opti¬ 
mal profile (I28|) . This minimization problem has a sim¬ 
ple mechanical analogy. Indeed, Eq. dSH) describes the 
motion of a classical particle in a constant gravity field 
(directed upward), and the solution with the jump of the 
first derivative corresponds to a bounce of the particle 
from the “floor” at h = M. Additional bounces (jumps 
in the hrst derivative) would cost more free energy, so 
they are not allowed. 

An extension of this theory to d dimensions is straight¬ 
forward. The short-time regime is a bit involved techni¬ 
cally, but the equilibrium (long-time) regime is simple. 
Here minimization of the free energy, constrained by the 
conservation law, yields a Gaussian distribution of M 
with the variance 


<Xd{L) 



(33) 


As one can see, d = 2 is the critical dimension [M3- 
The optimal surface at t = T exhibits a pyramid-like 
top. 

In conclusion, we have employed the MET to evalu¬ 
ate the probability that the stochastic Mullins-Herring 
interface with conserved noise first reaches a large given 
height at a specified time. We have also found the “opti¬ 
mal” time history of the interface conditional on the first 
passage. It would be interesting to apply the MET to 
find the range of the surface diffusion, that is to deter¬ 
mine the joint distribution of the (large) maximum and 
minimum of the interface height at a given time. The 
results of this work can be also applied to the stochas¬ 
tic Mullins-Herring equation with non-conserved noise, 
dfh = —d'^h + ^(x,t), with the same ^ as before. This 
equation describes, in a moving frarne, height fluctua¬ 
tions under particle deposition [^, [13 [I 4 I ■ By differenti¬ 
ating this equation with respect to x, one again arrives at 
(the one-dimensional version of) Eq. ([T]) but for the local 
slope of the interface s{x,t) = dxh{x,t). This mapping 
yields the probability that the interface first develops a 
(large) given local slope, and the optimal interface slope 
history conditional on this event. 

Einally, the exact solution of Eqs. o and m was pos¬ 
sible due to their linearity. Analogous MET equations 
for nonlinear stochastic equations, such as the Kardar- 
Parisi-Zhang equation (2^, are harder to solve. The sit¬ 
uation, however, is far from hopeless [^ . 
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Supplemental Material: Derivation of the MFT 
Equations 


The starting point of the derivation is the Langevin 
equation 

dth =-d^h + da;^{x,t), (34) 


where ^ is a delta-correlated Gaussian noise: 

= S{xi - X2)S{ti - t2). 


The problem is defined on the interval jxj < L/2 with 
periodic boundary conditions. We are interested in the 
probability of transition from the initial state h{x, t = 
0) = 0 to a given state h{x, t = T) at a specified time T. 
Because of the conservation law, h{x, T) must satisfy the 

condition T)dx = 0)dx = 0. Let us 

introduce a potential tp(x,t) so that dx'il’ = h{x,t). Now 
Eq. (m becomes 

dti’ = -dt'if-^{x,t), (35) 

and we can express 

^{x,t) =-dttp - (36) 


The Gaussian action is, therefore, 

pT 


S = 


/o 



nLl2 

dt 

/ dx 

J 

'-L/2 


rL/2 

dt 

/ dx 

J 

'-L/2 


(37) 


Being interested in large deviations, we need to minimize 
this action with respect to the trajectory 'ip(x,t). The 
variation of the action is 


6S = 


rL/2 


dt 


I-LI2 


dx{dtip + dx'4’){9tS'^ + d^Sip). (38) 


Let us introduce the momentum density field p{x,t), so 
that 


dxP = -dtip - (39) 


This follows 


dtxi^ = -dx{dxP + 


or 


dth = -dx{dlh + dxp) 


(40) 
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one of the Hamilton equations. Now we can rewrite the 
variation (1551) as follows: 

pT rL/2 

6s = — dt dx dxp{dtStjj + . 

Jo J-LI 2 

After several integrations by parts, we obtain the Euler- 
Lagrange equation 

dxtP = dip (41) 


which yields the second Hamilton equation of the MFT 
formalism: 

dtp = dtp. (42) 

The boundary terms in time and in space, resulting from 
the integrations by parts, all vanish because we specified 
fixed states at t = 0 and t = T, and because of the 
periodic boundary conditions in x. 


